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Abstract 

New, approximate, two-electron wavefunctions are introduced for 
the two-electron atoms (cations), which account remarkably well for 
the ground-state energies and the lowest- excited states (where avail- 
able). A new scheme of electronic configurations is also proposed for 
the multi-electron atoms. 

Following Bohr, [0| the electronic structure of the atoms is currently de- 
scribed by means of the central self-consistent field, based on the Hartree- 
Fock equations. H'[[| Although this mean- field doctrine has long served the 
interpretation of the atomic spectra, its central assumption of single-electron 
wavefunctions has never been given an a priori legitimacy. In fact, correlated 
ground-state wavefunctions have been worked out in great detail for the He 
atom and, generally, for the two-electron atom. [^]~ J/J In spite of the fact that 
these very elaborate techniques, which are variational in essence, produce im- 
pressively accurate numerical results for the ground-state energies, they look, 
however, rather arbitrary, and can not be extended easily to multi-electron 
atoms. We introduce in this paper new, approximate, two-electron wave- 
functions for the two-electron atom, which give the ground-state and the 
lowest-excited states energies with an accuracy which may be regarded as 
being remarkable for the simplicity of the approach. We propose also a new 



general scheme of electronic configurations, which may constitute a reason- 
able starting point for understanding the nature of the electronic structure 
of the multi-electron atoms. 

If not specified otherwise, we use throughout this paper the atomic units 
an — f^/me 2 = 0.53 A (Bohr radius) and e 2 /a# = 27.2 eV (twice the 
rydberg), which render the square of the electron charge e 2 = 1; in addition, 
we set the Planck's constant H = 1, so that the electron mass m — 1. We 
start with considering two electrons, denoted by 1 and 2, placed at ri and, 
respectively, r 2 in the Coulomb field of a nuclear charge Z. We neglect the 
center-of-mass corrections, the spin-orbit coupling, as well as all the other 
relativistic corrections; they can be treated perturbationally. Under these 
assumptions the hamiltonian of the two-electron atom is given by 

1 o 1 o Z Z 1 , . 

h = -pI + -pI +, 1 , 1 

2 2 ri r 2 | r x - r 2 | 

where pi i2 are the momenta of the two electrons. We shall treat first the 
electron affinity of hydrogen, i.e. the ground-state of the H~ anion (Z = 1), 
and thereafter we shall proceed to the atoms (cations) with Z > 2. 

If we put ri = — r 2 in (1) (for Z — 1) we get a one-particle hamiltonian 
for a particle of mass 1/2, moving in a Coulomb field of charge 3/2; its 
ground-state energy is Eq ^ = —9/16 = —15.3 eV, which amounts to an 
electron affinity = 1.7 eK.This value is close to the experimental value 
A = 0.76 eV.|| The zeroth order approximation Eq ^ must be corrected 
by higher-order contributions of the Coulomb repulsion V — 1/ | 1*1 — r 2 | 
between the two electrons. For a fixed r x we can see that this Coulomb 
repulsion describes, to the first order of approximation, harmonic oscillations 
of r 2 = (r 2 , # 2 , <£> 2 ) along 8 2 , around r 2 = — r^. The corresponding harmonic- 
oscillator potentials have, however, various curvatures, and we should define 
an average potential. Setting to describe the small oscillations of r 2 by the 
displacements p, i.e. putting r 2 = — rj + p, we may expand the Coulomb 
repulsion V as 



V=; r = ^- + ^jPi(coB0) + ... , (2) 

ri — r 2 I 2r 4r z 



where r =\ ri | = | r 2 |, Pi is the Legendre polynomial of the first rank, 
and 6 is the angle between p and ri. We define the average potential V 
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/ — -\ 1/2 

by V = I V ) , where the average in paranthesis is taken over all the 



orientations of p, so that, we get in the first approximation 

The first term in (3) has already been included in computing the zeroth order 
approximation; the second term in (3) is a harmonic-oscillator potential, 
which is completely determined by replacing r by its average value over the 
zeroth order ground-state r — 2. We obtain therefore the harmonic-oscillator 
frequency w = (24r 3 ) -1 ^ 2 = 1.96 eV, and the ground-state energy corrected 
by the zero-point oscillations E = + \w = —14.32 eV; this amounts 
to an electron affinity A = 0.72 eV which is in excellent agreement with the 
experimental value 0.76 eV. We remark that the oscillations of the second 
electron are along # 2 , s° that they correspond to only one degree of freedom. 

For Z > 2 we define r = min (r 1; r 2 ) and R = max (r 1? r 2 ) and separate 
the hamiltonian (1) as follows: 

H = H z + H Zl + V , (4) 

v =x^\-\ • (7) 

where Z\ — Z — 1 and p, P are the momenta of the electrons placed at 
r and, respectively, R. This separation corresponds to the classical con- 
cept of screening. The two hamiltonians H z and H Zl have a hydrogen-like 
energy spectrum; the ground-state wavefunction is given by ^o( r i) r 2) = 
Co^ r ioo( r )^ r ioo^(R')) where ^[qq Zi ^ are the corresponding ground-state wave- 
functions of the hydrogen-like atom (normalized to unity) and Co is a normal- 
ization constant. This wavefunction corresponds to an electron configuration 
which may be denoted by (100) x (100), and has the energy 

# = -Y-f = ~Z 2 + Z-l/2 . (8) 
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One can check easily that this energy is very close to the actual ground-state 
energy of the two-electron atoms (cations). For example, we have = 
—68 eV as compared with the experimental value Eq XP = —79 eV for the He 
atom (Z = 2), and E^ 0) = -176.8 eV, as compared with E™ p = -198 eV 
for the Li + cation. The remaining interaction V can be written as 

OO I OO I A I 

^ = E^(cosfi) = E E ^•^•*U^)^(©,<J>) , (9) 

1=1 1=1 m=—l 

where Pi are the Legendre polynomials, Yi m are the spherical harmonics, 
Q is the angle between r and R , and r = (r, #,</?), R = (R,©,$). This 
interaction can, in principle, be treated perturbationally. However, the main 
contribution to this interaction comes from those space regions where r ~ R, 
so that the perturbation series converges extremely slowly. Consequently, we 
shall account for its main effect in another way. 

Analyzing its expression given by (7) we see easily that V has a minimum 
value at r = R and Q = tc; this value is given by V = — l/2r. We shall 
assume, therefore, that the ground-state wavefunction given above is slightly 
distorted in such a way as to allow the second electron to take full advantage 
of this minimum of energy. This amounts to correcting E^ by V = — l/2fo, 
where r is the average value of the electron distances to the nucleus T\ = 
r 2 = r over the hydrogen-like configuration (100) x (100) given above. We 
give here the expression of this average radius 

r = Cl ■ 16(ZZ 1 ) 3 J dr 2 ■ r\ ■ {£ dn ■ r\ ■ e" 2Zri • e~ 2Z ^ + 

POO 

+ / dn-rf- e" 2Ziri • e~ 2Zr2 ) . (10) 

The ground-state energy is therefore 

E = E^ - l/2f , (11) 

and we remark that — l/2r is a classical correction, corresponding to the 
principle of minimization of the energy for the ground-state. In particular, we 
note that the first-order correction of V, as given by (9), within the quantum- 
mechanical perturbation theory vanishes for the ground-state. Within the 
present picture the ground-state of the two-electron atom looks as being 
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made of two electrons, one moving inside, and the other outside an imaginary 
sphere, and changing continuously their places through ri = — r 2 . Remark 
that further corrections like zero-point oscillations around r 2 = — r 1 done for 
the H~ anion would be inappropriate here since the second electron moves 
everywhere (independently) outside the r 2 = r\ sphere in the Coulomb field 
of charge Z\. The normalization constants Co, the average radii Tq and the 
ground-state energies E given by (11) are computed in Table 1 for Z = 2 
through 11, and compared with the experimental values of the ground-state 
energies Eq XP . One can se that the agreement is very good, the errors being 
less than 0.6%. 

Within the present approximation the full wavefunctions of the ground- 
state of the two-electron atom, including the spin degrees of freedom, may 
be written (up to a normalization constant) as 

' MSfa^fo) , n<r 2 , \ 



*o(ri,r 2 ) = [a(l)/3(2)-a(2)/3(l)]x 

for the spin-singlet state, and 
a(l)a(2) 



*SS(r 2 )*SS ) (n) , ri>r 2 



(12) 




*o(ri,r 2 ) = < 

/3(l)/3(2) J ^ _w iool r 2j^ioo 

(13) 

for the spin-triplet state, where a, f3 are the spin-up and spin-down wave- 
functions. The wavefunctions written above are antisymmetric, such as to 
satisfy the Pauli exclusion principle. Within the present approximation these 
wavefunctions correspond to the same energy, so that the para- and the ortho- 
atoms are degenerate. However, a closer inspection of (12), for example, tells 
us that the orbital of this wavefunction, though continuous at ri = r 2 has 
discontinuous derivatives at these points; and, similarly, the orbital wave- 
function given by (13) is discontinuous at = r 2 . The exact orbitals will be 
continuous, and will have continuous derivatives at = r 2 , which implies 
that the spin-triplet orbital will have a sudden variation on passing through 
zero at ri = r 2 . This indicates that the ortho-atom will have a higher energy, 
so that we may conclude that the ground-state corresponds to a spin-singlet 
(para-atom) . 

We remark also that the wavefunctions given by (12) and (13) are,essentially, 
correlated wavefunctions, or genuine two-electron wavefunctions. Indeed, one 
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may define, for example, 



100 1 

fSS'tn) , n > r 2 



and write the spin-singlet orbital of the ground-state as <^2(l)<£i(2); in spite of 
certain appearances, this is a correlated, two-electron wavefunction, since, for 
example, according to its definition (14), the waefunction (^2(1), as function 
of ri, depends on r 2 . 

We pass now to the lowest-excited states of the two-electron atoms (cations), 
with nuclear charge Z > 2. It is easily to check that, to the zeroth or- 
der approximation, the lowest-excited states correspond to the hydrogen-like 
configurations (100) x (200) and (100) x (21m), i.e. the wavefunctions are 

Moo( r )^2oo^(R) and M^M^imC^-)) according to the notations introduced 
above. The corresponding energy is 

Z 2 Z\ 5Z 2 Z 1 

E 1 = 1 - = + . 15 

2 8 8 4 8 v ; 

Indeed, one can check easily that the other alternative, corresponding to 
the configurations (200) x (100), (21m) x (100), has a much higher energy. 
One can also check straightforwardly that the energies given by (15) are 
already very close to the experimental values, as one can see in Table 2, for 
Z = 2 (He) and Z = 3 [Li + ); the experimental values for larger Z, where 
available, have no longer been included in Table 2 as they correspond clearly 
to much higher-excited states; for example, E exp = —119.6 eV for the Be 2+ 
cation, g while (15) gives E = -248.2 eV for Z = 4. In Table 2 there 
have also been included the normalization constants C\ for the configuration 
(100) x (21m), as well as the average radii f\ over this state for Z = 2 through 
11. We shall remark, first, that a correction of the type — l/2ri, as in the case 
of the ground-state energy, is not appropriate here, since the excited states 
are only stationary states, but they do not, of course, minimize the energy, as 
the ground-state does. One may check that including — l/2fi in the energy 
computed above would, indeed, result in worse numbers, i.e. there would be 
a larger discrepancy with respect to the experimental values. Secondly, we 
remark that although the two configurations (100) x (200) and (100) x (21m) 
are degenerate in this approximation, the former is not orthogonal to the 
ground-state, as the latter is. An orthogonalization procedure will push this 
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state toward higher energies, so that the lowest-excited state corresponds 
very closely to the (100) x (21m) configuration. Thirdly, we may also remark 
that an improved approximation to the lowest excited level may proceed by 
standard perturbation calculations applied to the interaction V given by (9), 
keeping in mind that we have to include those states that conserve the total 
angular momentum. In the case of (100) x (21m) this is L — 1, and the 
interaction V couples this state to (21m) x (100), (21m) x (200), (32m') x 
(21m), etc. As we have said, the perturbation series converges, however, 
very slowly, as the main contributions come from r ~ R. And finally, let us 
remark that the continuity of the two-particle orbitals at ri = and of their 
derivatives, brings an additional splitting in energy between the spin-singlet 
and spin-triplet states, as discussed above. 

Obviously, the above scheme of electronic configurations may be ex- 
tended to multi-electron atoms. For example, for an atom with N elec- 
trons, we may set up (to the zeroth order approximation) the configura- 
tion (nlm) x (n I m ) x (n I m ) x ... etc, where each paranthesis denotes a 
hydrogen-like orbital corresponding, respectively, to a nuclear charge Z, Z\ = 
Z — 1, Z 2 = Z — 2, etc, up to Zn-i — Z — (N — 1). The remaining inter- 
action of the type V given in (7) can then be minimized classically for the 
ground-state, according to the procedure described above. Similar electronic 
configurations may also be used as starting points for the excited states. Ob- 
viously, we have, in this picture, correlated, multi-electronic wavefunctions, 
though each electron moves in a central field, corresponding, however, to dif- 
ferent nuclear charges. One may infere that, within this picture, the atoms 
are more "rarefied", leading to increased transition probabilities, as com- 
pared with those corresponding to all electrons moving in a unique nuclear 
charge. Of course, the first thing to be done for checking the validity of such a 
scheme of electronic configurations would be that of trying to account for the 
periodicity of the chemical elements, i.e. the counterpart of the closed-shell 
assumption of the self-consistent field model of the atom. 
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Table 1 



z 


Co 


*"0 


-E (eV) 


-E^(eV) 


2 (He) 


0.795 


1.177 


79.554 


78.98 


3 (Li 1+ ) 


0.856 


0.650 


197.713 


198 


4 (Be 2+ ) 


0.889 


0.452 


370.087 


371.5 


5 (B 3+ ) 

6 (C 4 +) 


0.910 


0.347 


596.809 


599.3 


0.925 


0.282 


877.91 


881.6 


7 (iV 5+ ) 


0.935 


0.237 


1.213 • 10 3 


1.218 • 10 3 


8 (0 6+ ) 


0.943 


0.205 


1.603- 10 3 


1.609- 10 3 


9 (F 7 +) 


0.949 


0.180 


2.048 • 10 3 


2.055 • 10 3 


10 (iVe 8+ ) 


0.954 


0.161 


2.546 • 10 3 


2.556 • 10 3 


11 (iVa 9 +) 


0.958 


0.145 


3.099 • 10 3 


3.11 • 10 3 



Table 2 



z 




fi 


-E 1 (eV) 


-ET v {eV) 


2 (He) 


0.709 


2.907 


57.8 


59.16 


3 (Li 1+ ) 


0.712 


1.558 


136 


138.98 


4 (Be 2+ ) 


0.714 


1.092 


248.2 




5 (5 3 +) 

6 (C 4 +) 


0.715 


0.853 


394.4 




0.717 


0.708 


574.6 




7 (iV 5 +) 


0.718 


0.610 


788.8 




8 (0 6 +) 


0.718 


0.539 


1.037- 10 3 




9 (F 7 +) 
10 (A^e 8+ ) 


0.719 


0.485 


1.319 • 10 3 




0.719 


0.444 


1.635- 10 3 




11 (Na 9+ ) 


0.720 


0.410 


1.986- 10 3 
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Table captions 
Table 1 

The normalization constants C , the average radii r , the ground-state 
energies Eq as given by (11) and the experimental values Eq XP of the ground- 
state energies|| for the ground-state configuration (100) x (100) of the two- 
electron atoms (cations), as functions of the nuclear charge Z. 

Table 2 

The normalization constants C±, the average radii r 1; the energies E\ as 
given by (15) for the lowest-excited states corresponding to the configuration 
(100) x (21m), as well as the experimental values Ef xp for the lowest levels [§] 
of the two-electron atoms (cations), as functions of the nuclear charge Z; the 
experimental values for Z > 3 have not been included as they correspond to 
higher excited levels. 
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